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Velocity  Slip  and  Temperature  Jump  in  Hypersonic 

Aerothermodynamics 

Andrew  J.  Lofthouse*  Leonard  C.  Scalabrin*  and  Iain  D.  Boydt 

Department  of  Aerospace  Engineering,  University  of  Michigan,  Ann  Arbor,  4 8109 


Hypersonic  vehicles  experience  different  flow  regimes  during  flight  due  to  changes  in 
atmospheric  density.  Computational  Fluid  Dynamics  (CFD),  while  relatively  computa¬ 
tionally  inexpensive,  are  not  physically  accurate  in  areas  of  highly  non-equilibrium  flows. 
The  direct  simulation  Monte  Carlo  (DSMC)  method,  while  physically  accurate  for  all  flow 
regimes,  is  relatively  computationally  expensive.  A  breakdown  parameter  can  be  used 
to  determine  where  in  the  flow  domain  the  CFD  methods  are  valid.  The  current  study 
investigates  the  effect  of  continuum  breakdown  on  surface  aerothermodynamic  properties 
(pressure,  shear  stress  and  heat  transfer  rate)  of  a  cylinder  in  Mach  10  and  Mach  25  flows 
of  argon  gas  for  several  different  flow  regimes,  from  the  continuum  to  a  rarefied  gas.  Sev¬ 
eral  different  velocity  slip  and  temperature  jump  boundary  conditions  are  examined  for 
use  with  the  CFD  method.  CFD  and  DSMC  solutions  are  obtained  at  each  condition. 
Total  drag  and  peak  heat  transfer  rate  predictions  by  CFD  remains  within  about  6%  of  the 
DSMC  predictions  for  all  regimes  considered,  with  the  generalized  slip  condition  proposed 
by  Gokgen  giving  the  best  results. 


Nomenclature 


A  Constant 

d  Molecular  diameter 

k  Boltzmann  constant 

L  Characteristic  length 

m  Mass  of  one  molecule 

q  Heat  transfer  rate 

T  Translational  Temperature 

Tj  Temperature  jump 

Tw  Wall  temperature 

Velocity  in  x-direction 
Slip  velocity 
p  Pressure 

Symbols 

X  Mean  free  path 

p  Coefficient  of  viscosity 

u  VHS  temperature  exponent 

p  Mass  density 

a  Momentum  accommodation  coefficient 

r  Shear  stress  parallel  to  wall 
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THE  VIEWS  EXPRESSED  IN  THIS  ARTICLE  ARE 
THOSE  OF  THE  AUTHOR  AND  DO  NOT  REFLECT 
THE  OFFICIAL  POLICY  OR  POSITION  OF  THE 
UNITED  STATES  AIR  FORCE,  DEPARTMENT  OF 
DEFENSE,  OR  THE  U.S.  GOVERNMENT. 


I.  Introduction 


The  design  of  hypersonic  vehicles  requires  accurate  prediction  of  the  surface  properties  while  in  flight. 
These  quantities  are  typically  the  heat  flux,  pressure  and  shear  stress.  These  variables  govern  not  only  the 
aerodynamic  performance  of  the  vehicle,  but  also  determine  the  selection  and  sizing  of  the  thermal  protection 
system  (TPS),  which  protects  the  vehicle  from  the  entry  environment.  During  it’s  trajectory  through  an 
atmosphere,  a  hypersonic  vehicle  will  experience  vastly  different  flow  regimes  because  the  atmosphere’s 
density  varies  as  a  function  of  altitude.  Reproduction  of  these  varied  flow  conditions  in  ground-based 
laboratory  facilities  is  both  expensive  and  technically  challenging.  Hence,  there  is  an  extremely  important 
role  for  computational  models  in  the  development  of  hypersonic  vehicles. 

At  low  altitudes,  the  atmospheric  density  is  relatively  high,  and  flows  around  hypersonic  vehicles  should  be 
simulated  using  traditional  Computational  Fluid  Dynamics  (CFD)  by  solving  either  the  Euler  or  preferably 
the  Navier-Stokes  (NS)  equations.  This  is  the  continuum  regime  characterized  by  very  large  Reynolds 
numbers  and  very  low  Knudsen  numbers,  where  the  Knudsen  number  is  defined  in  Eq.  1 . 


„  A  1 
Kn  =  —  oc  — 
L  pL 


a) 


At  very  high  altitudes,  at  the  edge  of  the  atmosphere,  the  density  is  low  such  that  there  are  very  few 
collisions  between  the  molecules  and  atoms  in  the  flow  around  the  vehicle.  This  is  the  rarefied  flow  regime 
and  can  be  computed  using  the  direct  simulation  Monte  Carlo  (DSMC)  method.1  Generally  speaking,  CFD 
methods  for  solving  the  NS  equations  are  about  an  order  of  magnitude  faster  than  the  DSMC  method. 
However,  the  lack  of  collisions  makes  the  physics  of  the  NS  equations  invalid  in  rarefied  regimes,  which  is 
characterized  by  a  large  Knudsen  number.  On  a  blunt  body,  a  high-density  fore-body  flow  can  create  a 
rarefied  flow  in  the  wake  of  the  vehicle.  In  principle,  the  DSMC  method  can  be  applied  to  any  dilute  gas 
flow,  but  becomes  prohibitively  expensive  for  Knudsen  numbers  less  than  0.001.  Thus,  it  is  attractive  to  find 
ways  to  increase  the  validity  of  CFD  methods  beyond  the  continuum  regime. 

One  way  to  improve  CFD  modeling  in  the  transition  regime,  that  is,  for  higher  density  flows  beyond 
the  continuum  regime  and  before  the  free  molecular  regime,  is  by  replacing  the  typical  no-slip  boundary 
conditions  with  slip  velocity  and  temperature  jump  boundary  conditions.  The  addition  of  slip  boundary 
conditions  will  not,  however,  eliminate  all  source  of  errors  when  using  continuum  methods  for  flows  with 
large  amounts  of  non-equilibrium.  For  instance,  in  more  rarefied  flows,  the  pressure  tensor  ceases  to  be 
isotropic,  a  key  assumption  in  the  derivation  of  the  continuum  equations. 

The  areas  of  the  flow  where  the  continuum  hypothesis  breaks  down  (or  equivalently,  where  the  flow  is 
no  longer  in  local  thermodynamic  equilibrium),  can  be  quantified  by  the  use  of  a  continuum  breakdown 
parameter.  While  there  have  been  several  breakdown  parameters  presented  in  the  literature,2  5  the  one 
that  appears  to  be  most  appropriate  for  hypersonic  compressible  flows  is  the  gradient-length  local  (GLL) 
Knudsen  number6,7 


KncLL  — 


A_ 

Q 


dQ 

dl 


(2) 


where  the  derivative  is  taken  in  the  direction  of  the  maximum  gradient,  and  Q  is  some  quantity  of  interest 
such  as  density,  pressure,  temperature  or  velocity  magnitude.  When  calculating  KncLL  based  on  velocity 
magnitude,  the  gradient  is  normalized  by  the  maximum  of  the  local  velocity  magnitude  and  the  local  speed 
of  sound.  It  is  generally  assumed  that  continuum  breakdown  occurs  whenever  KncLL  (based  on  the  CFD 
solution)  is  greater  than  0.05. 

Previous  work8  sought  to  quantify  the  errors  in  CFD  predictions  of  surface  properties  of  a  hypersonic 
blunt  body  by  comparing  those  solutions  with  DSMC  solutions  for  a  variety  of  flows  in  the  transitional  regime. 
The  maximum  KncLL  was  also  computed  in  each  case  to  determine  the  effect  of  continuum  breakdown  on 
the  surface  properties.  The  cases  considered  were  limited  to  a  Mach  10  flow  of  argon.  The  previous  CFD 
solutions  were  obtained  with  no-slip  boundary  conditions  imposed  at  the  wall.  As  would  be  expected  when 
only  considering  a  simple  gas,  the  results  indicated  that  the  inclusion  of  velocity  slip  and  temperature  jump 
boundary  conditions  would  improve  the  agreement  between  the  DSMC  and  CFD  conditions. 

The  current  work  builds  on  these  results  with  the  main  objective  being  to  quantify  the  improvement 
in  agreement  between  the  DSMC  solutions  and  the  CFD  solutions  with  the  inclusion  of  slip  boundary 
conditions.  Additional  cases  of  a  Mach  25  flow  are  also  computed  to  determine  any  effects  of  continuum 
breakdown  at  higher  velocities  on  the  surface  properties. 
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This  paper  will  first  briefly  describe  the  simulation  procedures  used,  including  some  background  of  the 
computational  models.  Several  slip  boundary  conditions  are  then  evaluated.  It  will  then  discuss  general 
flow  field  and  surface  property  results  predicted  by  the  models  in  several  different  flow  regimes,  from  the 
continuum  to  a  rarefied  gas.  Finally,  some  conclusions  are  presented. 

II.  Background  and  Simulation  Procedure 

This  investigation  continues  previous  work  that  considered  a  hypersonic  flow  of  argon  over  a  two- 
dimensional,  12-inch  diameter  cylinder,  as  shown  in  Figure  1.  Two  difference  free  stream  velocities  are 
used  that  correspond  to  Mach  10  and  Mach  25  flows  (UQ 0  =  2624.1  m/s  and  6585  m/s,  respectively).  The 
wall  temperature  is  held  constant  at  500  K  (Mach  10)  and  1500  K  (Mach  25).  The  free  stream  density  of  the 
flow  is  varied  such  that  several  different  regimes  are  considered,  from  the  continuum  through  the  transitional 
to  the  rarefied  regime,  as  shown  in  Table  1.  Knudsen  numbers  are  calculated  based  on  freestream  conditions 
and  the  cylinder  diameter,  using  the  hard-sphere  model  for  the  mean  free  path  calculation.  Surface  and  flow 
field  properties  for  this  flow  are  presented  from  two  different  computational  approaches. 


Figure  1.  Geometry  definition. 


First,  CFD  results  are  obtained  through  solu¬ 
tion  of  the  Navier-Stokes  equations.  The  CFD  re¬ 
sults  are  obtained  using  the  Michigan  Aerothermo- 
dynamic  Navier-Stokes  (LeMANS)  code,  developed 
at  the  University  of  Michigan  for  the  simulation  of 
hypersonic  reacting  flow- fields.9, 10  LeMANS  is  a 
general  2D/3D,  parallel,  unstructured,  finite- volume 
CFD  code  capable  of  simulating  gases  in  thermal 
and  chemical  non-equilibrium.  All  LeMANS  solu¬ 
tions  are  generated  assuming  an  isothermal  wall  at 
500  K  (Mach  10)  and  1500  K  (Mach  25).  Different 
boundary  conditions,  including  no-slip  and  slip  ve¬ 
locity  and  temperature  jump  are  enforced.  These 
will  be  discussed  in  more  detail  below.  In  each  case, 
a  grid  independence  study  was  conducted  to  determine  the  final  mesh  resolution  used. 

Second,  DSMC  results  are  provided  from  the  MONACO  code11  for  the  same  flow  conditions.  MONACO 
is  a  general  2D/3D,  parallel,  unstructured  mesh  DSMC  code  that  has  been  applied  to  many  hypersonic, 
rarefied  flows.12  All  MONACO  solutions  are  generated  using  a  fixed  wall  temperature  at  500  K  and  1500  K. 
Bird’s  variable  hard  sphere  model  is  used.1  In  general,  the  mesh  used  for  the  final  solution  for  each  case  is 
adapted  from  previous  solutions  such  that  each  cell  size  is  on  the  order  of  a  mean  free  path.  The  exceptions 
are  the  Kn  =  0.002,  where  the  cell  size  is  approximately  four  times  the  mean  free  path,  and  the  Kn  =  0.01, 
where  the  cell  sizes  near  the  surface  are  on  the  order  of  two  mean  free  paths  in  size.  In  these  cases  the  subcell 
method  is  used  to  select  particles  for  collisions1  to  ensure  physical  accuracy. 

The  variable  hard  sphere  (VHS)  parameters  used  in  the  DSMC  computations  are  equivalent  to  those  used 
in  the  previous  work  (the  temperature  exponent,  cj,  is  0.734  with  a  reference  diameter  of  3.595  x  10“ 10  m  at 


Table  1.  Flow  regimes  considered. 


Krioo’ 

Mass  Density 
(kg/m3) 

Number  Density 
(particles/m3) 

0.002 

0.01 

0.05 

0.25 

1.408  x  10"4 
2.818  x  lO"5 
5.636  x  10“6 
1.127  x  10-6 

2.124  x  1021 
4.247  x  1020 
8.494  x  1019 
1.699  x  1019 

*  Based  on  hard-sphere  mean  free  path. 
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a  reference  temperature  of  1000  K).  The  viscosity  and  thermal  conductivity  of  pure  argon  in  the  CFD  code 
are  computed  to  exactly  correspond  to  the  VHS  model,  according  to  Eqs.  3  and  4. 13 

(3) 

1 5 \J  7rinkTre  / 

/ire/  =  2nd2ref  (5  -  2a ;)  (7  -  2u)  [  } 

Some  computational  details  are  given  in  Table  2.  It  should  be  noted  here  that  both  MONACO  and 
LeMANS  are  research  codes.  Additional  improvements  to  DSMC  codes  have  been  proposed14  that  may 
reduce  the  computational  expense  of  the  DSMC  calculations.  The  LeMANS  code  is  very  new  and  has  not 
been  extensively  optimized  for  performance.  Also,  the  CFD  cases  cited  in  Table  2  correspond  to  the  Type 
2  slip  condition  (discussed  below)  which  had  the  slowest  convergence  rate.  Additional  work  will  most  likely 
yield  much  better  computational  efficiency. 

Table  2.  DSMC  and  CFD  computational  details. 


Mach  10 


Knoo 

DSMC 

CFD 

Cells 

Particles 

Time  Steps 

CPUs 

Cells 

Iterations 

CPUs 

0.002 

394,250 

287.0x10® 

1,050,000 

128 

80,000 

60,000 

8 

0.01 

68,497 

26.8x10® 

430,000 

32 

40,000 

55,000 

8 

0.05 

18,818 

7.1x10® 

430,000 

32 

40,000 

87,000 

8 

0.25 

24,452 

8.2x10® 

290,000 

32 

40,000 

90,000 

8 

Mach  25 


Knoo 

DSMC 

CFD 

Cells 

Particles 

Time  Steps 

CPUs 

Cells 

Iterations 

CPUs 

0.002 

121,063 

142.0x10® 

1,400,000 

128 

80,000 

42,000 

8 

0.01 

87,364 

32.6x10® 

460,000 

32 

40,000 

50,000 

8 

0.05 

33,268 

18.7x10® 

450,000 

32 

40,000 

65,000 

8 

0.25 

14,006 

9.7x10® 

300,000 

32 

40,000 

67,000 

8 

II. A.  Slip  Boundary  Conditions 

This  section  discusses  several  different  types  of  slip  boundary  condition  and  their  implementation  in  the 
LeMANS  code. 

The  slip  boundary  condition  was  first  derived  by  Maxwell  for  a  flat  plat.  The  expression  he  derived  for 
the  slip  velocity  is  shown  in  Eq.  5. 15 


Va  =  A 


3  jj_0T 

4  pT  dx 


(5) 


For  an  isothermal  wall,  the  temperature  gradient  can  be  neglected,  giving  the  boundary  condition  in  its 
simplest  form  in  Eq.  6. 

/  O  \ 

(6) 


The  mean  free  path,  A,  is  calculated  from  typical  gas  flow  properties  as  shown  in  Eq.  7. 


16 


2/i  =  M  /  7T 

pc  p  V  2  RT 


(7) 
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The  boundary  condition  for  temperature  jump  is  similar,  seen  in  Eq.  8. 17 


To -Tw  = 


2  —  a  27  ^  dT 

a  (7  -I-  1)  Pr  dy 


0 


(8) 


The  thermal  accommodation  coefficient,  a,  can  be  assumed  to  be  the  same  as  the  momentum  accommodation 
coefficient,  a.  In  this  case,  the  second  and  third  terms  on  the  right  hand  side  can  be  combined  and  an 
equivalent  mean  free  path  adopted  as 


At  = 


4  K 
(7  +  1)  pcCv 


2  /C  /  7T 

(7+l)paV2flT 


to  give  the  simpler  form  shown  in  Eq.  10. 


To  -Tw  = 


2  —  a  dT 
- At 

a  dy 


(9) 


(10) 


The  simple  Maxwell  slip  conditions  given  by  Eqs.  6  and  10  are  implemented  in  LeMANS  and  referenced 
below  as  the  Type  1  slip  condition. 

Throughout  many  years  of  research,  the  general  form  of  the  velocity  slip  has  remained  the  same,  while 
the  difficulty  is  in  determining  the  correct  value  for  the  constant,  A,  and  the  accommodation  coefficients,  a 
and  a.  Many  times  these  coefficients  are  determined  empirically,  while  other  times  they  are  calculated  from 
relations  based  on  kinetic  theory.18  The  most  detailed  boundary  conditions  require  knowledge  a  priori  of  a 
good  non-equilibrium  solution  for  the  flow  being  considered.  In  effect,  any  change  in  A  simply  changes  the 
value  of  the  accommodation  coefficient,  which  depends  on  the  physical  characteristics  of  the  wall  itself.19 
The  current  investigation  seeks  to  find  a  good  boundary  condition  that  does  not  require  a  priori  knowledge 
of  any  non-equilibrium  solution,  which  would  eliminate  any  advantage  in  using  CFD  over  DSMC.  Also, 
both  the  DSMC  and  CFD  solutions  are  calculated  assuming  a  fully  diffuse  wall,  which  corresponds  to  an 
accommodation  coefficient  of  1. 

The  previous,  simple  slip  boundary  conditions  were  derived  for  small  Knudsen  numbers.  Gokcen16  showed 
that,  for  large  Knudsen  numbers,  these  simple  slip  boundary  conditions  converge  to  different  values  than 
those  predicted  by  free  molecular  flow.  He  then  proposed  the  general  slip  boundary  conditions  given  in  Eq. 


11. 


ota  ( Q>\  Uu»)  —  2Aa 


(11) 


Here,  a  is  either  velocity  or  temperature,  and  Aa  is  given  by  either  Eq.  7  or  9.  For  small  Knudsen  numbers, 
this  expression  reduces  to  Eqs.  6  and  10. 17  This  slip  condition  is  referenced  in  the  results  below  as  the  Type 
2  slip  condition. 

As  mentioned  previously,  the  original  Maxwell  slip  boundary  condition  was  derived  for  a  flat  plate. 
Lockerby  et  al.15  showed  that  for  a  curved  wall  the  proper  boundary  condition  is  given  by  Eq.  12.  This  is 
known  as  the  extended  slip  boundary  condition. 


2-a\  ( du *  dun\  3  p  dT 
a  )  \  dn  dx  )  4  pT  dx 


(12) 


This  slip  condition  is  not  implemented  in  LeMANS  at  this  time  and  is  not  considered  further  here. 

The  Navier-Stokes  equations  fail  in  at  least  two  different  areas  as  the  flow  becomes  more  rarefied.  The  first 
area,  an  assumption  of  no-slip  flow,  is  corrected  through  the  use  of  the  slip  boundary  conditions.  However, 
the  N-S  equations  also  assume  that  the  shear  stress  varies  linearly  with  the  velocity  gradient  at  the  wall. 
Cercignani  showed  the  form  that  the  velocity  profile  in  the  Knudsen  layer  takes.20 

Here  there  are  two  choices  one  can  make  concerning  the  boundary  conditions;  use  the  actual  slip  velocity 
in  which  case  the  flow  field  away  from  the  wall  will  not  be  accurate,  or  use  a  fictitious  slip  velocity.  The 
fictitious  slip  will  not  accurately  predict  the  flow  near  the  wall  (and  hence  may  adversely  affect  the  accuracy 
of  drag  and  heat  transfer  rates)  but  will  be  more  accurate  in  the  region  away  from  the  wall. 

A  correction  to  the  velocity  gradient  is  evaluated  by  Lockerby  et  al.20  They  proposed  using  a  wall- 
function  type  of  boundary  condition  for  the  velocity  where  the  viscosity  in  the  Knudsen  layer  is  modified 
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a s  fi  =  ^  1  where  ^  is  given  in  Eq.  13,  derived  from  a  curve- fit  to  Cercignani’s  Knudsen  layer  velocity 
profile. 


This  modification  of  the  viscosity  in  the  Knudsen  layer  is  also  used  in  connection  with  a  simple  slip  boundary 
condition  as  given  by  Eq.  6  with  A  =  y/n/2.  This  approach  is  expected  to  allow  the  CFD  method  to 
accurately  model  the  velocity  profile  in  the  Knudsen  layer. 

This  new  wall  function  has  only  been  evaluated  for  isothermal  flow  conditions  that  are  usually  encountered 
in  micro-flows.  Although  there  will  be  some  change  to  the  heat  transfer  rate  coefficient  (based  on  a  constant 
Prandtl  number),  there  might  be  other  changes  necessary  to  give  the  correct  temperature  profile  in  the 
Knudsen  layer  for  the  non-isothermal,  hypersonic  flows  considered  here.  Nevertheless,  it  is  instructive  to 
investigate  the  possible  improvement  this  wall  function  might  afford.  This  slip  condition  is  implemented  in 
LeMANS  and  is  referenced  as  the  Type  3  slip  condition. 

The  traditional  no-slip  boundary  conditions  are  referenced  below  as  the  Type  0  boundary  conditions. 

II. B.  Slip  Quantities  and  DSMC 

The  Boltzmann  Equation  inherently  includes  slip  boundary  conditions;  therefore,  the  DSMC  method  includes 
them  as  well.  In  MONACO,  the  slip  velocity  and  temperature  jump  are  calculated  based  on  the  particles  that 
strike  the  surface,  according  to  Eqs.  14  and  15,  which  are  the  relations  used  in  Bird’s  DSMC  implementation, 
DS2V. 21,22 


y  =EWW 
3  EWK) 

E  (m/Vn)  (t u 2  +v2  +  w2)-Y:  (m/Vn)V? 
3*E(1/V„) 


(14) 

(15) 


Comparison  of  MONACO  and  DS2V  results  show  that  the  velocity  slip  and  temperature  jump  results 
(using  Eqs.  14  and  15)  are  nearly  identical.  This  is  not  surprising  considering  that  their  implementations 
are  identical.  However,  it  is  useful  to  compare  the  DSMC  slip  velocity  and  temperature  jump  to  those 
obtained  by  the  CFD  slip  conditions  previously  mentioned,  in  order  to  evaluate  their  relative  merit.  It 
appears,  however,  that  the  quantities  computed  according  to  Eqs.  14  and  15  may  not  be  the  best  choice 
when  comparing  solutions.  Figure  2  compares  the  CFD  and  DSMC  data  at  the  wall  boundary  for  the 
Mach  10,  Kn  =  0.01  case.  The  DSMC  data  is  given  by  Eqs.  14  and  15  while  the  CFD  data  is  given  by 
the  slip  velocity  and  temperature  jump  according  the  equations  in  the  previous  section,  depending  on  the 
case.  Figure  3  compares  the  same  case,  but  the  data  is  extracted  from  the  field  solution  5  x  10-7  m  from 
the  surface  of  the  cylinder,  which  is  well  within  the  Knudsen  layer.  Although  Figure  2  suggests  that  the 
CFD  solution  predicts  a  higher  velocity  slip  and  temperature  jump  (and  the  temperature  jump  doesn’t  even 
follow  the  same  general  trend),  Figure  3  suggests  that  the  agreement  is  actually  much  closer,  especially  for 
the  Type  2  boundary  condition.  Therefore,  it  appears  that  data  extracted  from  the  field  data  might  be  more 
appropriate  to  compare  than  the  actual  boundary  data. 


III.  Results 

This  investigation  seeks  to  quantify  the  differences  in  the  surface  properties  of  a  cylinder  in  a  hypersonic 
flow  by  comparing  CFD  and  DSMC  predictions  for  total  drag  and  peak  heat  transfer  .  Additionally,  the 
overall  flow  features,  such  as  the  temperature  field,  are  compared,  as  well  as  the  surface  distributions  of 
pressure,  shear  stress  and  heat  transfer. 

Tables  3  and  4  compare  the  total  drag  and  the  peak  heat  transfer  rates  predicted  by  both  computational 
methods.  Here,  the  peak  heat  transfer  rate  is  obtained  by  averaging  over  the  surface  area  within  one  degree  of 
the  stagnation  point.  For  CFD,  these  quantities  are  calculated  for  each  of  the  different  boundary  conditions 
implemented,  with  Type  0  being  the  no-slip  condition. 

It  is  clear  that  the  slip  boundary  conditions  improve  the  agreement  between  the  two  simulation  methods. 
While  the  difference  in  total  drag  is  as  high  as  23.6%  and  32.0%  for  the  no-slip  condition,  the  Type  2  slip 
boundary  condition  shows  the  best  agreement,  with  the  difference  exceeding  5%  only  for  the  Kn  =  0.25  case. 
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-  MONACO  -  MONACO 

- CFD(O)  - CFD(O) 


(a)  Velocity  Slip  (b)  Temperature  Jump 

Figure  2.  Comparison  of  flow  properties  at  the  wall  for  CFD  and  DSMC  using  slip  boundary  conditions  (CFD)  and 
Eqs.  14  and  15  (DSMC).  The  DSMC  properties  are  significantly  less  than  all  the  CFD  solutions. 


-  DSMC 

- CFD  (0) 


-  DSMC 

- CFD  (0) 


(a)  Velocity  Magnitude  (b)  Temperature 

Figure  3.  Comparison  of  flow  properties  for  CFD  and  DSMC  by  extracting  data  from  the  field  properties  5  x  10“ 7 m 
from  the  wall.  Agreement  between  the  two  methods  is  better. 
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Similar  trends  are  noted  for  the  peak  heat  transfer  rate,  with  a  nearly  40%  difference  for  the  no-slip  CFD 
solution,  while  the  Type  2  slip  solution  remains  near  1%  or  less  for  most  of  the  cases  and  only  exceeds  5% 
for  the  Kn  =  0.25  case  at  Mach  25. 

There  is  also  a  slight  improvement  in  drag  and  heat  flux  agreement  between  the  Type  1  and  Type  3 
boundary  conditions,  where  the  only  difference  in  the  two  boundary  conditions  is  the  inclusion  in  the  Type 
3  of  the  viscosity  correction  function  within  the  Knudsen  layer.  This  improvement  is  not  as  great  as  that 
achieved  with  the  Type  2  boundary  condition,  but  it  might  hint  that  the  Type  2  condition  could  be  improved 
with  the  addition  of  the  viscosity  correction  in  the  Knudsen  layer. 

Since  the  Type  2  condition  shows  the  best  agreement  with  the  DSMC  solutions,  the  CFD  temperature 
field  results  shown  below  are  taken  from  the  Type  2  CFD  case. 

Table  3.  Total  drag. 


Mach  10 


Kn^ 

Drag/Length  [N/m]  (%  Difference) 

DSMC 

CFD  (0) 

CFD  (1) 

CFD  (2) 

CFD  (3) 

0.002 

187.6 

187.5  (-0.1%) 

187.4  (-0.1%) 

187.4  (-0.1%) 

187.6  (0.0%) 

0.01 

40.02 

40.32  (0.8%) 

40.19  (0.4%) 

40.16  (0.4%) 

40.27  (0.6%) 

0.05 

8.900 

9.416  (5.8%) 

9.122  (2.5%) 

8.863  (-0.4%) 

9.078  (2.0%) 

0.25 

2.092 

2.585  (23.6%) 

2.301  (10.0%) 

1.982  (-5.2%) 

2.253  (7.7%) 

Mach  25  • 


Knoo 

Drag/Length  [N/m]  (%  Difference) 

DSMC 

CFD  (0) 

CFD  (1) 

CFD  (2) 

CFD  (3) 

0.002 

0.01 

0.05 

0.25 

1,171 

250.8 

56.85 

13.35 

1,176  (0.4%) 
255.3  (1.8%) 
61.17  (7.6%) 
17.61  (32.0%) 

1,176  (0.4%) 
254.6  (1.5%) 
59.38  (4.5%) 
15.73  (17.8%) 

1,176  (0.4%) 
254.4  (1.4%) 
57.04  (0.3%) 
13.02  (-2.4%) 

1,177  (0.5%) 
255.1  (1.7%) 
58.90  (3.6%) 
15.14  (13.4%) 

The  breakdown  parameter  is  calculated  using  both  the  CFD  and  the  DSMC  solutions  using  Equation 
2  with  the  derivative  being  taken  in  the  direction  of  the  steepest  gradient.  For  this  case  of  a  cylinder  in 
a  hypersonic  flow  of  a  simple  gas,  the  only  causes  of  breakdown  to  the  continuum  hypothesis  expected  are 
in  regions  of  high  gradients  (such  as  the  shock  and  boundary  layer)  and  regions  of  rarefaction  (such  as  in 
the  wake).  The  amount  of  continuum  breakdown  is  also  expected  to  increase  as  the  gas  flow  becomes  more 
rarefied.  These  trends  are  confirmed  by  the  results  that  are  shown  below.  In  general,  the  flow  experiences 
continuum  breakdown  in  three  areas;  across  the  bow  shock,  in  the  boundary  layer  and  in  the  wake  region. 
The  flow  in  the  shock  and  boundary  layer  regions  experiences  very  steep  gradients  in  flow  properties,  while 
the  wake  region  is  more  rarefied,  thus  leading  to  the  breakdown  of  the  continuum  hypothesis. 


In  the  results  that  follow,  the  surface  properties  are  presented  in  terms  of  non-dimensionalized  coefficients 
with  the  following  definitions: 


^  P-Poo 

P  \PooUl 

(16) 

°F  ~  ipM 

(17) 

Ch  =  I  p  u> 

2  Poo  ^  oo 

(18) 

where  p  is  the  pressure,  r  is  the  shear  stress,  q  is  the  heat  transfer  rate,  poc 

is  the  free  stream  pressure,  p^ 

is  the  free  stream  density  and  is  the  free  stream  velocity.  The  surface  properties  in  each  case  are  plotted 
as  a  function  of  the  angle  around  the  cylinder,  with  the  stagnation  point  being  located  at  an  angle  of  zero 
(see  Figure  1). 
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Table  4.  Peak  heat  transfer  rate. 


Mach  10 


Krioo 

Peak  Heating  [kW/m2] 

(%  Difference) 

DSMC 

CFD  (0) 

CFD  (1) 

CFD  (2) 

CFD  (3) 

0.002 

89.80 

89.85  (0.0%) 

89.14  (-0.7%) 

89.00  (-0.9%) 

88.54  (-1.4%) 

0.01 

39.13 

40.22  (2.8%) 

39.49  (0.9%) 

39.20  (0.2%) 

38.89  (-0.6%) 

0.05 

15.92 

18.08  (13.6%) 

17.25  (8.4%) 

15.79  (-0.8%) 

16.81  (5.6%) 

0.25 

5.926 

7.851  (32.5%) 

7.061  (19.1%) 

6.184  (4.3%) 

6.640  (12.1%) 

Mach  25 


Kn^ 

Peak  Heating  [kW/m2]  (%  Difference) 

DSMC 

CFD  (0) 

CFD  (1) 

CFD  (2) 

CFD  (3) 

0.002 

0.01 

0.05 

0.25 

1,746 

749.6 

309.9 

105.7 

1,763  (0.9%) 
791.0  (5.5%) 
357.9  (13.6%) 
147.8  (39.9%) 

1,750  (0.2%) 
777.4  (3.7%) 
341.0  (10.4%) 
133.6  (19.1%) 

1,746  (0.0%) 
768.4  (2.5%) 
294.2  (-4.8%) 
111.7  (5.7%) 

1,739  (-0.4%) 
765.3  (2.1%) 
331.2  (7.2%) 
123.7  (17.0%) 

Along  with  the  surface  properties,  the  maximum  KncLL  at  the  surface  (based  on  the  CFD  solution  with 
Type  2  slip  conditions)  is  also  plotted  in  each  case. 


IV.  Kn  =  0.002 

At  a  Knudsen  number  of  0.002,  the  flow  is  within  the  continuum,  no-slip  regime.  Nevertheless,  there  is 
still  evidence  of  continuum  breakdown  in  the  shock,  along  the  cylinder  surface  in  the  boundary  layer  and  in 
the  wake  as  seen  in  Figure  4.  Interestingly,  DSMC  predicts  a  larger  degree  of  breakdown  than  does  CFD. 
The  increase  in  velocity  from  Mach  10  to  Mach  25  also  slightly  increases  the  amount  of  breakdown  predicted 
in  the  shock  and  wake  regions. 


Figure  4.  Kn^  =  0.002  KriGLL  for  Mach  10  and  Mach  25.  Light  gray  regions  correspond  to  KncLL  <  0.05,  dark  gray 
regions  correspond  to  KnGLL  <  0.10  and  black  regions  correspond  to  KnGLL  >  0.10.  Continuum  breakdown  is  predicted 
in  the  shock  region,  in  a  thin  boundary  layer  along  the  surface  and  in  a  region  of  flow  expansion  around  the  top  of  the 
cylinder. 
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The  temperature  field  predicted  by  CFD  and  DSMC,  as  seen  in  Figure  5,  is  in  excellent  agreement,  as 
expected.  The  Mach  25  case  shows  an  increase  in  the  maximum  temperature  as  well  as  a  slight  increase 
in  the  shock  thickness;  however,  for  the  presently  considered  simple  gas,  there  are  no  other  effects.  The 
surface  properties  shown  in  Figures  6-8  show  very  good  agreement  between  the  two  methods  for  all  types 
of  CFD  boundary  conditions  considered,  and  for  both  Mach  10  and  Mach  25.  The  only  location  along  the 
surface  where  KnQLL  is  greater  than  0.05  is  in  the  wake  region.  The  total  drag  and  peak  heat  transfer  rates 
predicted  by  the  two  methods  differ  by  less  than  1%,  with  only  one  exception  (that  being  the  Type  3  CFD 
solution  at  1.4%).  These  results  are  well  within  the  error  range  expected  for  a  flow  for  which  both  methods 
are  valid. 


(a)  Mach  10 


ICFDl 


(b)  Mach  25 


Figure  5.  Krioo  =  0.002  temperature  field,  for  Mach  10  and  Mach  25. 


£ 

E 


<t>  (deg) 


MONACO 


(a)  Mach  10 


(b)  Mach  25 


Figure  6.  Knoo  =  0.002  surface  pressure  coefficient  and  KnGLL  for  Mach  10  and  Mach  25. 
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£ 

£ 


(a)  Mach  10  (b)  Mach  25 

Figure  7.  Krioo  =  0.002  surface  friction  coefficient  and  KriGLL  for  Mach  10  and  Mach  25. 
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Figure  8.  Kn*  =  0.002  surface  heating  coefficient  and  Kiigll  for  Mach  10  and  Mach  25. 
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Maximum  Knou  Maximum  Knolx 


V.  Kn  =  0.01 


(a)  Mach  10  (b)  Mach  25 

Figure  9.  Kn^,  =  0.01  KnGLL  for  Mach  10  and  Mach  25.  Light  gray  regions  correspond  to  KnGLL  <  0.05,  dark  gray 
regions  correspond  to  KnGLL  <  0.10  and  black  regions  correspond  to  KnGLL  >  0.10.  There  are  larger  regions  where 
continuum  breakdown  becomes  a  concern,  especially  in  the  wake,  the  shock  and  the  boundary  layer.  The  Mach  25  case 
shows  larger  regions  than  the  Mach  10  case. 


A  Knudsen  number  of  0.01  is  considered  to  be  near  the  limit  of  the  continuum,  no-slip  regime.  Here  there 
is  increased  evidence  of  continuum  breakdown,  as  shown  in  Figure  9.  However,  comparison  of  general  flow 
field  features  shows  that  the  two  numerical  solutions  (CFD  and  DSMC)  do  not  differ  much.  The  temperature 
field,  seen  in  Figure  10,  is  very  similar,  with  a  few  exceptions  in  the  shock  thickness  and  the  wake  (for  Mach 
25),  where  the  continuum  hypothesis  is  expected  to  break  down  first.  The  shock  standoff  distance  predicted 
by  both  methods  is  the  same,  as  is  the  maximum  temperature  behind  the  shock,  for  both  Mach  10  and  Mach 
25.  The  thicker  shock  does  not  seem  to  have  an  effect  on  the  surface  properties,  as  is  shown  below. 


(a)  Mach  10  (b)  Mach  25 


Figure  10.  Krioo  =  0.01  temperature  field  for  Mach  10  and  Mach  25. 


The  amount  of  breakdown  in  the  shock  does  not  necessarily  carry  over  to  the  surface  properties.  Again, 
the  amount  of  continuum  breakdown  near  the  surface  does  not  seem  to  significantly  affect  the  surface 
properties,  shown  in  Figures  11  -  13.  The  no-slip  CFD  case  (Type  0)  shows  some  divergence  from  the  others 
in  the  wake.  The  best  agreement  is  achieved  with  the  Type  2  and  Type  3  slip  conditions.  Shear  stress 
starts  to  diverge,  even  with  the  slip  conditions,  at  Mach  25.  This  particular  case  seems  to  be  somewhat  of 
an  aberration,  especially  when  considering  the  good  agreement  that  will  be  shown  below  for  the  Kn  =  0.05 


12  of  23 


American  Institute  of  Aeronautics  and  Astronautics  Paper  2007-0208 


case.  The  heating  rates  also  agree  well  at  Mach  10,  but  show  some  differences  a  the  higher  velocities  of  Mach 
25.  The  total  drag  due  to  pressure  and  viscous  effects  predicted  by  CFD  is  still  within  1%  of  that  predicted 
by  DSMC,  for  Mach  10,  and  within  2%  for  Mach  25,  as  shown  in  Table  3.  The  peak  heating  also  differs  by 
less  than  2%  for  the  Mach  10  case  when  using  slip  boundary  conditions,  and  less  than  3%  for  the  best  slip 
conditions  (Type  2  and  Type  3)  for  the  Mach  25  case.  These  results  indicate  that  the  surface  properties 
are  not  strongly  affected  by  the  continuum  breakdown  in  the  shock,  and  the  slip  conditions  improve  the 
agreement  quite  satisfactorily. 


MONACO 


MONACO 


(a)  Mach  10  (b)  Mach  25 

Figure  11.  Kn^  =  0.01  surface  pressure  coefficient  and  KnGLL  for  Mach  10  and  Mach  25. 


MONACO 


MONACO 


(a)  Mach  10  (b)  Mach  25 

Figure  12.  Kn^  ==  0.01  surface  friction  coefficient  and  KncLL  for  Mach  10  and  Mach  25. 
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MONACO 


MONACO 


(a)  Mach  10 


(b)  Mach  25 


Figure  13.  Kn^ 


=  0.01  surface  heating  coefficient  and  KncLL  for  Mach  10 


and  Mach  25. 
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Maximum  Kn. 


VI.  Kn  =  0.05 


ICFDl 


fCFBI 


(a)  Mach  10 


(b)  Mach  25 


Figure  14.  Kn^o  =  0.05  KnGLL  for  Mach  10  and  Mach  25.  Light  gray  regions  correspond  to  KnGLL  <  0.05,  dark  gray 
regions  correspond  to  KnGLL  <  0.10  and  black  regions  correspond  to  KnGLL  >  0.10.  Continuum  breakdown  is  predicted 
in  much  larger  regions  than  previous  cases.  The  DSMC  method  predicts  larger  regions  where  continuum  breakdown 
exceeds  the  critical  value  of  0.05,  with  an  even  larger  region  predicted  for  the  Mach  25  case,  cylinder. 


At  a  Knudsen  number  of  0.05,  the  flow  is  well  within  the  slip  regime.  The  flow  demonstrates  breakdown 
in  a  larger  area  of  the  flow,  in  each  of  the  three  regions  (shock,  boundary  layer  and  wake),  as  seen  in  Figure 
14.  Again,  there  are  larger  areas  in  the  Mach  25  case  where  the  breakdown  parameter  exceeds  the  critical 
value  of  0.05;  in  fact,  nearly  the  entire  domain  is  predicted  to  be  beyond  the  region  of  applicability  of  the 
continuum  equations.  The  differences  between  the  CFD  and  DSMC  temperature  fields  shown  in  Figure  15 
are  more  pronounced  than  the  lower  Knudsen  number  cases,  with  even  more  differences  for  the  higher  Mach 
number  case.  The  DSMC  shock  is  much  thicker  than  the  CFD  shock,  although  the  shock  stand-off  distance 
and  peak  temperatures  are  still  very  nearly  equal.  In  the  wake,  DSMC  predicts  a  higher  temperature  than 
CFD. 


fCFOl 


fCFDl 


(a)  Mach  10 


(b)  Mach  25 


Figure  15.  Kiioo  =  0.05  temperature  field  for  Mach  10  and  Mach  25. 


The  surface  pressure  predicted  by  both  methods  is  still  in  excellent  agreement,  as  seen  in  Figure  16. 
However,  the  shear  stress  is  higher  for  most  of  the  CFD  cases,  especially  in  the  wake,  with  the  surprising 
exception  of  the  Type  2  case,  which  agrees  almost  exactly  with  the  DSMC  case  (see  Figure  17).  The  total 
drag  for  all  of  the  slip  cases  are  within  5%  of  the  DSMC  case,  with  the  Type  2  case  being  almost  identical. 
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The  heat  transfer  rate,  shown  in  Figure  18,  also  shows  surprisingly  good  agreement  between  the  Type  2  case 
and  DSMC.  While  the  no-slip  case  predicts  a  peak  heat  transfer  rate  about  13%  higher  than  DSMC,  the  slip 
cases  show  better  agreement,  with  the  Type  2  case  within  1%  for  Mach  10  and  less  than  5%  for  Mach  25. 


MONACO 


MONACO 


(a)  Mach  10  (b)  Mach  25 

Figure  16.  Kn^  =  0.05  surface  pressure  coefficient  and  Kiigll  for  Mach  10  and  Mach  25. 
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Figure  17.  Kn*  =  0.05  surface  friction  coefficient  and  Kiigll  for  Mach  10  and  Mach  25. 


The  fortuitous  agreement  between  the  Type  2  CFD  case  and  the  DSMC  case  raises  the  question  of  how 
good  the  agreement  is  between  flow  properties,  such  as  temperature  and  velocity,  further  from  the  surface. 
Recall  that  using  the  actual  slip  velocity  and  temperature  jump  as  boundary  conditions  in  a  continuum 
method  can  improve  wall  properties  such  as  shear  stress  and  heat  transfer,  but  may  introduce  additional 
errors  in  the  predicted  flow  properties  further  from  the  wall.  To  investigate  this  possibility,  the  velocity  and 
temperature  (along  with  the  breakdown  parameter)  are  plotted  along  a  line  normal  to  the  surface  at  an  angle 
of  90°,  and  is  shown  in  Figures  19  and  20.  Both  plots  show  that  there  is  very  good  agreement  between  the 
CFD  slip  cases  and  the  DSMC  solution  near  the  wall  and  extending  out  into  the  flow  field  for  the  Mach  10 
case.  For  Mach  25,  there  is  some  disagreement  in  the  temperature  at  the  wall,  although  it  recovers  further 
out,  while  the  agreement  in  velocity  magnitude  predicted  is  not  as  good.  Therefore,  there  doesn’t  seem  to 
be  any  concern  about  significant  errors  in  flow  properties  predicted,  at  least  in  the  case  of  the  Mach  10  flow. 
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Figure  18.  Krioo  =  0.05  surface  heating  coefficient  and  KriGLL  for  Mach  10  and  Mach  25. 
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Figure  19.  Kn,*  =  0.05  velocity  magnitude  and  KnGLL  along  a  line  normal  to  body  surface  at  4>  =  90°.  The  region  of 
continuum  breakdown  spreads  further  from  the  surface.  A  velocity  slip  of  about  200  m/s  is  predicted  by  DSMC. 
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Figure  20.  Kn^  =  0.05  velocity  magnitude  and  KncLL  along  a  line  normal  to  body  surface  at  4>  =  90°.  The  region  of 
continuum  breakdown  spreads  further  from  the  surface.  A  velocity  slip  of  about  200  m/s  is  predicted  by  DSMC. 


VII.  Kn  =  0.25 
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Figure  21.  Kn  =  0.25  KncLL  for  Mach  10  and  Mach  25.  Light  gray  regions  correspond  to  KncLL  <  0.05,  dark  gray 
regions  correspond  to  KnGLL  <  0.10  and  black  regions  correspond  to  KnGLL  >  0.10.  Continuum  breakdown  is  predicted 
in  almost  the  entire  domain  surrounding  the  cylinder. 


For  Knoo  =  0.25,  the  flow  would  be  considered  outside  of  the  slip  regime  and  into  the  transition  regime. 
Here,  even  the  addition  of  slip  boundary  conditions  is  not  expected  to  help  the  continuum  CFD  method’s 
predictive  capabilities  very  much.  Indeed,  the  plots  of  the  breakdown  parameter  in  Figure  21  indicate  that 
there  are  significant  non-equilibrium  effects  across  almost  all  of  the  flow  domain.  The  temperature  fields 
predicted  by  both  methods,  seen  in  Figure  22,  shows  some  major  differences.  Although  the  peak  temperatures 
behind  the  shock  appear  to  be  about  the  same,  the  shock  stand-off  distances  are  significantly  different,  as  is 
the  thermal  boundary  layer  thickness. 

The  surface  pressure,  plotted  in  Figure  23,  predicted  by  both  methods  is  also  no  longer  in  agreement, 
even  for  the  CFD  cases  implementing  slip  conditions.  The  DSMC  pressure  is  less  than  the  CFD  pressure 
near  the  fore- body  (most  likely  due  to  the  thinner  shock  predicted  by  CFD,  which  allows  the  CFD  solution  to 
approach  the  ideal  jump  condition  more  closely  than  the  DSMC  solution),  but  the  agreement  does  improve 
in  the  wake.  The  pressure  tensor  at  this  Knudsen  number  is  also  most  likely  non-isotropic,  which  would 
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Figure  22.  Kn  =  0.25  temperature  field  for  Mach  10  and  Mach  25. 


introduce  errors  into  the  CFD  solution.  The  shear  stress,  shown  in  Figure  24,  shows  the  same  general  trend 
as  in  previous  cases  in  that  both  methods  agree  near  the  stagnation  region,  but  the  results  diverge  as  the 
flow  accelerates  around  the  cylinder.  Here,  the  Type  2  CFD  solution  predicts  a  lower  shear  stress,  but  a 
higher  pressure,  thus  the  overall  drag  predictions  are  somewhat  close  (differing  by  about  5%  for  Mach  10 
and  less  than  3%  for  Mach  25).  The  heat  transfer  rate,  as  seen  in  Figure  25,  also  follows  trends  similar  to 
the  previous  cases  in  that  the  DSMC  heat  transfer  rate  is  lower  than  the  CFD  rate  along  the  entire  surface 
for  all  but  the  Type  2  case,  which  gets  fairly  close  in  some  areas.  The  peak  heat  transfer  rates  differ  by  as 
much  as  almost  40%  for  the  no-slip  case  at  Mach  25,  but  the  Type  2  case  shows  the  best  agreement  with 
less  than  5%  difference  at  Mach  10  and  less  than  6%  difference  at  Mach  25. 
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Figure  23.  Krioo  =  0.25  surface  pressure  coefficient  and  Kiigll  for  Mach  10  and  Mach  25. 


The  velocity,  temperature  and  breakdown  parameter  plotted  along  the  line  normal  to  the  surface  at  an 
angle  of  90°  are  shown  in  Figures  26  and  27.  In  this  case,  there  are  significant  differences  between  the  Type 
2  case  and  DSMC,  while  the  Type  1  and  Type  3  cases  show  better  agreement.  It  should  be  noted  that  the 
cells  near  the  surface  of  the  DSMC  cases  might  not  be  small  enough  to  capture  the  gradients  in  the  Knudsen 
layer. 
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Figure  24.  Kn^  =  0.25  surface  friction  coefficient  and  KnGLL  for  Mach  10  and  Mach  25. 
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Figure  25.  Knw  =  0.25  surface  heating  coefficient  and  KncLL  for  Mach  10  and  Mach  25. 
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Figure  26. 


Knoo  =  0.25  velocity  magnitude  and  Kiigll  along  a  line  normal  to  body  surface  at  4> 


90°. 
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Figure  27.  Kn^  =  0.25  temperature  and  KnGLL  along  a  line  normal  to  body  surface  at  $ 


90°. 
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VIII.  Conclusion 


The  primary  goal  of  this  study  was  to  quantify  the  difference  between  DSMC  and  CFD  simulations, 
when  using  slip  boundary  conditions,  to  determine  the  effects  of  different  levels  of  continuum  breakdown  on 
surface  properties  under  hypersonic  flow  conditions. 

Comparison  of  CFD  and  DSMC  results  for  identical  flow  conditions  showed  that  the  surface  properties 
of  pressure,  shear  stress  and  heat  transfer  rates  were  very  similar  for  the  lower  Knudsen  number  flows  where 
the  continuum  hypothesis  is  valid,  as  expected,  while  the  results  diverged  in  the  higher  Knudsen  number 
cases.  The  surface  pressure  was  least  affected  by  continuum  breakdown,  as  quantified  by  the  gradient-length 
local  Knudsen  number,  among  those  properties  investigated,  and  seemed  to  be  affected  only  by  breakdown 
in  the  shock  region  at  the  highest  Knudsen  number  flow.  The  shear  stress  was  most  influenced  by  continuum 
breakdown  and  was  affected  primarily  by  continuum  breakdown  in  the  wake.  The  addition  of  slip  velocity 
and  temperature  jump  boundary  conditions  greatly  improved  the  agreement  at  higher  Knudsen  numbers. 
Several  different  types  of  slip  boundary  conditions  were  examined,  and  the  best  agreement  appears  to  be 
obtained  when  using  the  generalized  slip  conditions  proposed  by  Gokgen.  With  these  boundary  conditions, 
the  difference  in  surface  properties  predicted  by  CFD  and  DSMC  increased  from  less  than  1%  at  Krioo  =  0.002 
to  around  5%  at  Kn^  =  0.25. 

For  the  case  of  a  simple  gas,  the  higher  velocities  associated  with  a  Mach  25  flow  did  not  seem  to  increase 
the  difference  between  the  CFD  and  DSMC  predictions.  Although  the  extent  of  the  region  where  the 
continuum  breakdown  parameter  exceeded  the  critical  value  of  0.05  was  larger  at  the  higher  Mach  number, 
the  predicted  surface  properties  with  the  slip  boundary  conditions  stilled  remained  well  under  5%  for  all  but 
the  Kn^  =  0.25  case,  where  the  peak  heat  transfer  rates  and  total  drag  predictions  were  within  6%. 

Future  Work 

The  present  case  was  limited  to  a  simple  gas.  Gases  comprised  of  diatomic  molecules,  such  as  nitrogen, 
are  expected  to  present  additional  difficulties,  such  as  vibrational  relaxation  and  dissociation  rates,  especially 
for  the  higher  temperatures  encountered  at  Mach  25. 

Additionally,  different  geometries,  such  as  sharp  leading  edges  with  attached  shocks,  are  expected  to 
present  more  challenges  in  getting  CFD  predictions  to  agree  wtih  those  from  DSMC  due  to  the  highly 
non-equilibrium  behavior  near  the  leading  edge. 
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